Duality, Magnetic space group and their applications to quantum phases and phase 
transitions on bipartite lattices in several experimental systems 
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By using a dual vortex method, we study phases such as superfluid, solids, supersolids and 
quantum phase transitions in a unified scheme in extended boson Hubbard models at and slightly 
away from half filling on bipartite optical lattices such as honeycomb and square lattice. We also 
map out its global phase diagram at T = of chemical potential versus the ratio of kinetic energy 
over the interaction. We stress the importance of the self-consistence condition on the saddle point 
structure of the dual gauge fields in the translational symmetry breaking insulating sides, especially 
in the charge density wave side. We find that in the translational symmetry breaking side, different 
kinds of supersolids are generic possible states slightly away from half filling. We propose a new kind 
of supersolid: valence bond supersolid ( VB-SS). In this VB-SS, the density fluctuation at any site is 
very large indicating its superfluid nature, but the boson kinetic energies on bonds between two sites 
are given and break the lattice translational symmetries indicating its valence bound nature. We 
show that the quantum phase transitions from solids to supersolids driven by a chemical potential 
are in the same universality class as that from a Mott insulator to a superfluid, therefore have 
exact exponents z — 2, v = 1/2, 77 = with a logarithmic correction. Comparisons with previous 
quantum Monte-Carlo (QMC) simulations on a square lattice are made. Implications on possible 
future QMC simulations in both bipartite lattices are given. All these phases and phase transitions 
can be potentially realized in ultra-cold atoms loaded on optical bipartite lattices. Then we apply 
our results to investigate the reentrant "superfluid" in a narrow region of coverages in the second 
layer of 4 He adsorbed on graphite and the low temperature phase diagram of Hydrogen physisorbed 
on Krypton- preplated graphite ( f/2/Kr/graphite ) near half filling. We suggest that 4 He and H2 
lattice supersolids maybe responsible for the experimental signals in the two systems. Finally, we 
suggest Cooper supersolid is repressible for the phase diagram of La2- x Ba x Cu04 near x = 1/8. 



I. INTRODUCTION. 

The Boson Hubbard model with various kinds of in- 
teractions, on all kinds of lattices and at different filling 
factors is described by the following Hamiltonian 1 : 
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where rii — b\bi is the boson density, t is the nearest 
neighbor hopping which is determined by the depth of 
the trapping potential at the prefered adsorption sites. 
U, Vi , V2 are onsite, nearest neighbor (nn) and next near- 
est neighbor (nnn) interactions respectively, the • • • may 
include further neighbor interactions and possible ring- 
exchange interactions. 

In the hard-core limit U — > 00, due to the exact map- 
ping between the boson operator and the spin s = 1/2 op- 
erator: b\ = Sf, hi = S^,rii = 5? + 1/2, the boson model 
EqnJT] can be mapped to a "generalized" anisotropic 
5 = 1/2 quantum Heisenberg model in an external mag- 
netic neld±&^: 



H = -2t E (SfSJ + SfS*) 
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where h = fi — 2Vi — 2V2 for a square lattice and the 
• ■ • may include further neighbor interactions and ring- 
exchange interactions. 

The model Eqn[T]with only the onsite interaction was 
first studied in Ref.— . It was found that there is a sec- 
ond order superconductor to insulator transition at filling 
factor / = 1. The effects of long range Coulomb inter- 
actions on the transition was studied ir&£. The duality 
transformation to the dual vortex picture at / = n was 
performed ir4 and is briefly reviewed in the appendix. 
Here, we only outline its essence to facilitate the duality 
transformation in the more general cases / = p/q to be 
discussed in the next several paragraphs. In the direct 
boson picture, a vortex is a singularity in the boson wave- 
function, so a boson wavefunction acquires a 2-7T phase 
when in encircles a vortex. In the dual vortex picture, 
a boson is a singularity in the vortex wavefunction, so a 
vortex wavefunction acquires a 2-7T phase when in encir- 
cles a boson. After performing the boson-vortex duality 
transformation, the authors in^ obtained a dual theory 
of Eqn[T]in term of the interacting vortices ip hopping on 
the dual lattice subject to a fluctuating dual " magnetic 
field" . The average strength of the dual " magnetic field 
" A M through a dual plaquette is equal to the boson den- 
sity / = n, because 27rn is equivalent to 0, so the average 
value can be simply taken to be zero. It is important 
to stress that the average density of bosons is the same 
in both the SF and the Mott insulating side, namely, it 
takes the integer n on both sides, so the average strength 
of the dual magnetic field can be taken as zero on both 
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side, then the fluctuations in the dual gauge field reflects 
the fluctuations of the boson density. In the continuum 
limit, just from the gauge invariance, the final effective 
theory in terms of the vortex order parameter tp is just 
the scaler electrodynamics of the dual vortex tp coupled 
to the fluctuating gauge field described by Eqn lA2l 
In the superfluid state < tp >= 0, the gauge field is gap- 
less, while in the insulating state < tp >^ 0, the gauge 
field acquires a mass due to the Higgs mechanism. This 
duality transformation has been confirmed by Quantum 
Monte- Carlo (QMC ) in the first reference ini. 

Obviously, it is important to study other commensu- 
rate filling factors / = p/q ( p,q are relative prime num- 
bers ), because in addition to the superfluid phase, many 
other phases such as charge density wave or valence bond 
solid can only appear when q > 2. So novel phase transi- 
tions such as SF to CDW or SF to VBS can be realized 
only in such cases. Recently, the dual vortex method 
( DVM ) developed in^ was greatly expanded to study 
EqnQ]on a square lattice just at such generic filling fac- 
tors ir£. It turns out that the DVM at q > 2 is much 
more involved than the simplest 5=1 case due to the cru- 
cial involvement of the non-commutative projective space 
group at q > 2. The general procedures are the follow- 
ing. After performing the similar charge-vortex duality 
transformation as ins, the authors in 5 obtained a dual 
theory of EqnQ] in term of the interacting vortices tp a 
hopping on the dual lattice subject to a fluctuating dual 
" magnetic field". The average strength of the dual " 
magnetic field " through a dual plaquette is equal to the 
boson density f — p/q. This is similar to the Hofstadtcr 
problem of electrons moving in a crystal lattice in the 
presence of a magnetic field^i The projective represen- 
tation of the space group (PSG) dictates that there are 
at least q-fold degenerate minima in the mean field en- 
ergy spectrum ( when q = 1, the PSG just reduces to the 
usual commutative space group ). Near the supercon- 
ductor to the insulator transition, the most important 
tp a fluctuations will be near these q minima which can 
be labeled as tpi, I = 0, 1, • • • , q — 1 which forms a q di- 
mensional representation of the PSG. In the continuum 
limit, the final effective theory in terms of these q order 
parameters should be invariant under this PSG. In the 
superfluid state < tpi >= for every /, while in the insu- 
lating state < tpi > for at least one /. In the insulating 
state, there must exist some kinds of charge density wave 
(CDW) or valence bond solid ( VBS) states which may 
be stabilized by longer range interactions or possible ring 
exchange interactions included in EqnQ] The CDW or 
VBS order parameter was constructed to be the most 
general bilinear and gauge invariant combinations of the 

In a recent unpublished preprint^, the author stud- 
ied all the possible phases and phase transitions in the 
EBHM slightly away from half filling q = 2 on bipar- 
tite lattices such as honeycomb and square lattice. It 
was found ir*i£ that the dual vortex method at d = 2 
can achieve many important results which are very diffi- 



cult to achieve from the direct boson picture. Although 
square lattice has been studied by various analytic and 
numerical methods before, the boson Hubbard model on 
a honeycomb lattice was not studied in any details by 
both analytic and Quantum Monte-Carlo (QMC) meth- 
ods before the preprint^. It is not a Bravais lattice, so 
may show some different properties than those in square 
lattice. Experimentally, the honeycomb lattice could also 
be easily realized in ultra-cold atomic experiments to be 
discussed in section II. The honeycomb lattice is also the 
relevant lattice for adatom adsorption on substrates to 
be discussed in section VIII. I found that a supersolid ( 
SS ) state exists only away from half filling. A super- 
solid in EqnQ] is a state with both superfluid and solid 
order. A supersolid in Eqnf2] is defined as the simulta- 
neous orderings of ferromagnet in the XY component ( 
namely, < hi >^ ) and CDW in the Z component. Re- 
cently, by using the torsional oscillator measurement, a 
PSU's group lead by Chan observed a marked 1 ~ 2% su- 
perfluid component even in bulk solid 4 He at ~ 0.2K S . 
If this experimental observation indicates the existence 
of 4 He supersolid remains controversial 9 . However, it 
was established by spin wave expansion 12 and quantum 
Monte-carlo (QMC) 13 ' 14 ' 15 simulations that a supersolid 
state could exist in an extended boson Hubbard model 
(EBHM) with suitable lattice structures, filling factors, 
interaction ranges and strengths. 

It was also explicitly pointed out iniS that the DVM 
developed in 5 holds only in the superfluid ( SF ) and the 
valence bond solid ( VBS ) side where the saddle point of 
the dual gauge field can be taken as uniform, however, it 
fails in the charge density wave ( CDW ) side where the 
saddle point of the dual gauge field can not be taken as 
uniform anymore. So not only the fluctuations, but also 
the average values of the dual gauge fields are different 
on both sides. This is in sharp contrast to the super- 
fluid to the Mott transition at and near f = n described 
by Eqns lA2l and IA4I in the appendix. As explained be- 
low Eqn lA21 the average density of bosons is the same 
in both the SF and the Mott insulating side, namely, it 
takes the integer n on both sides, so the average strength 
of the dual magnetic field can be taken as a uniform zero 
on both sides, of course, the fluctuations of the gauge 
field are completely different on both sides. However, in 
the CDW side which breaks the translational symmetry, 
special care is needed to choose a correct saddle point of 
the dual gauge field in the CDW side to make the the- 
ory self-consistent, so a different action is needed in the 
CDW sideiS. This paper is an expanded version of the 
unpublished preprint^. In this expanded version, (1) by 
pushing the DVM to slightly away from commensurate 
filling factors and (2) also extending the DVM explicitly 
to the lattice symmetry breaking CDW side by choos- 
ing the corresponding self-consistent saddle points of the 
dual gauge field, I will map out the global phase diagram 
of the EBHM on bipartite optical lattices such as hon- 
eycome and square lattice at and near half filling and 
also investigate superfluid, solid, especially supersolids 
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and quantum phase transitions in the phase diagram in 
a unified scheme. 

The DVM is a magnetic space group ( MSG 
symmetry-based approach which, in principle, can be 
used to classify all the possible phases and phase tran- 
sitions after choosing correct saddle points for the dual 
gauge fields. But the question if a particular phase will 
appear or not as a ground state depends on the spe- 
cific values of all the possible parameters in the EBHM 
in Eqn[TJ so it can only be addressed by a microscopic 
approach such as Quantum Monte-Carlo (QMC) simula- 
tions. The DVM can guide the QMC to search for partic- 
ular phases and phase transitions in a specific model. Fi- 
nite size scalings in QMC in a specific microscopic model 
can be used to confirm the phases and the universality 
classes of phase transitions discovered by the DVM. The 
two methods are complementary to each other and both 
are needed to completely understand phases and phase 
transitions in Eqn[TJ 

The rest of the paper is organized as following, In sec- 
tion II, I will explicitly derive the generators of the mag- 
netic space group ( also called projective space group ) 
for q = 2 in a honeycomb lattice, construct the effective 
actions which are invariant under this MSG, also write 
down both the charge density wave and the valence bond 
order parameters to characterize the symmetry breaking 
sides in the insulating regimes. In section III and IV, I 
will describe the phases and phase transitions driven by 
the competition of kinetic energy and repulsive potential 
energy in Fig. 4 ( represented by the parameter r in 
Fig. 2 and Fig. 3 ) along the horizontal axis at the com- 
mensurate fillings and the phases and phase transitions 
driven by the chemical potential /i along the vertical axis 
slightly away from the commensurate fillings in Ising and 
easy-plane limit respectively. We also stress that satis- 
fying the self-consist condition for the dual gauge field 
in the Ising limit is absolutely necessary to achieve cor- 
rect answers in the CDW side. In section V, we study 
the EBHM in a square lattice and also discuss the still 
disputed so called deconfined quantum critical point in 
the easy-plane limit. In section VI, we compare our re- 
sults achieved by the DVM with some available QMC 
results and make implications on possible future QMC 
simulations in both lattices. Then we study the applica- 
tion of our results on 3 different experimental systems: 
ultra-cold atoms with long range interactions on opti- 
cal lattices in both square and honeycomb lattice in sec- 
tion VII, adatom adsorption on different substrates such 
as possible reentrant supersolids in the second layer of 
4 iJe adsorbed on graphite and Hydrogen adsorbed on 
Krypton-preplated graphite in honeycomb lattice in sec- 
tion VIII, possible cooper pair supersolid in high temper- 
ature superconductor La 2 - x Ba x CuOi in square lattice 
in section IX. finally, we reach conclusions in section X. 
In the appendix, we review the boson-vortex duality at 
integer filling / = n in which I will stress the role of the 
dual gauge field and its connection to the more general 
case of / = p/q developed 

ii&io 

and the main text of this 



paper. 



II. MAGNETIC SPACE GROUP, EFFECTIVE 
ACTION AND ORDER PARAMETERS IN THE 
DUAL VORTEX PICTURE. 

In this section, we will extend the DVM in 5 to study 
the EBHM EqnQ] in honeycomb lattice at and slightly 
away from q = 2. Honeycomb lattice ( solid line in 
Fig.l ) is not a Bravais lattice, so may show some dif- 
ferent properties than a square lattice. The dual lat- 
tice of the honeycomb lattice is a triangular lattice ( 
dashed line in Fig.l ). Two basis vectors of a primi- 
tive unit cell of the triangular lattice can be chosen as 
di = x, <?2 = — \x + ¥f-y, 3d = (?i + a 2 as shown in Fig.l. 
The reciprocal lattice of a triangular is also a triangular 
lattice and spanned by two basis vectors k = kibi + k 2 b 2 
with bi = x + -^,h = -^y satisfying 6, • Sj = 5 lj . The 
point group of a triangular lattice is Cq v ~ D§ which 
contains 12 elements. The two generators can be cho- 
sen as Cq = i?7r/3 , h ■ The space group also includes 
the two translation operators T\ and T 2 along a\ and 
a 2 directions respectively. The 3 translation operators 
T\, T 2 , Td, the rotation operator R-k/3, the 3 reflection op- 
erators Ii 7 I 2 ,Id, the two rotation operators around the 
direct lattice points A and B: ^,13,^3 of the MSG 
are worked out in the following. 

In the Landau gauge A — (0,Hx), the mean field 
Hamiltonian for the vortices hopping in a triangular lat- 
tice in the presence of / flux quanta per triangle in the 
tight-binding limit is^ : 

H 0v = — t'S^[\x + a\ X x\ + h.c. 

X 

+ \x + a 2 > e l27r2fai < x\ + h.c. 

+ \x + a d > e ^2/(a 1+ i/2) <S \ + h c ^ (3) 

where x = a\ai + a 2 a 2 denotes lattice points of the tri- 
angular lattice. Note that the total vortex Hamiltonian 
H v = Hq v + V where V is the interaction between vor- 
tices. Because V does not contain any Aharonov-Bohm 
(AB) phase factor from the non-trivial A^ background, 
so the magnetic space group is completely determined by 
the vortex kinetic term Hq v . V always commutes with 
the generators of the MSG given by Eqns. 14151101 So 
when constructing the representation of the MSG, we 
can ignore the V term without losing any generality. 

The Ti,T 2 ,Td and the rotation operator 11^/3 of the 
PSG are worked out ir£. Here we listed them in slightly 
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different notations: 

Ti = ^\S+S 1 > e l27r2/a2 <x\ 

X 

t 2 = \x + a2 >< x\ 

X 

^/3 = ]Te l27r/(a? ~ 2aia2) |ai-a2,ai ><ai,a 2 |(4) 

.r 

It can be shown that they all commute with H . How- 
ever, they do not commute with each other T\T 2 — 
ujT d ,TxT 2 = u 2 T 2 T 1 . 

After performing some algebras, we found three of the 
6 reflection operators in the point group C§ v : 

h = if^e~ j27r/a 2|ai -a 2 ,-a 2 >< o 1 ,a 2 | 
h = K^2e~ l2Trfa2l \ - ax, aa - a x >< a x ,a 2 | 

X 

I d = /v^e- j27r/2aiQ2 |a 2 ,a 1 >< ai ,a 2 | (5) 

X 

where K is the complex conjugate operator which make 
Ix,I%,Is non-unitary operators. Note that in contrast 
to the reflection operators in square lattice, the phase 
factors in Eqn[5]are crucial to ensure they commute with 
the Hamiltonian Eqnj3] 

The eigenvalue equation Hip(k) = E(k)ip(k) leads to 
the Harper's equation in the triangle lattice : 

(e-iki + e -<(fei+fc2+27r/(2*-i))^_ i ( A . i) k2) 

+ {e tkl + e i ^ 1+fe2+27r/ ^- 1 »)^+i(^i, k 2 ) 
+2cos(-k 2 + 2nfl)Mki,k 2 )=E(k)^ l (kx,k 2 ) (6) 

where I = 0, 1, • • ■ , q — 1. 

In the following, we focus on q = 2 case where there is 
only one band E(k) — — 2i(cosfci+cosA; 2 — cos(fci +fc 2 )). 
Obviously, E(kx,k 2 ) = E(-k u -k 2 ) = E(k 2 ,kx). There 
are two minima at k± = ±(7r/3, tt/3). Let's label the 
two eigenmodes at the two minima as By using the 
expressions of the rotation and reflection operators in 
Eqns. UEl we find the two fields transform as: 

Tx,T 2 : V±^e™/ 3 V±; T d : V± -> -e^/ 3 ^ 
-Rtt/3 : "0± -* ^=f! 7 " : ^± a = 1, 2, d (7) 

where the transformations under T a , Rn/3 were already 
derived in 5 . Note that R^/?, plays the same role as the 
Z 2 exchange symmetry between tp + and ^>_. 

The quadratic terms of the effective action is simply 
the scalar electrodynamics as in the square lattice case: 

£o=EH- *AM 2 + r|Vg 2 + \Fl v (8) 



where — \e^ v \d v A\ is the strength of the non- 
compact gauge field A^. 

It is easy to show that there are only 2 independent 
quartic invariants under the above transformations: 

A 1 = hM 4 + hM 4 , A 2 = |^+| 2 kM 2 (9) 

The direct honeycomb lattice is a non-Bravais lattice 
which contains two lattice points A and B per direct unit 
cell, it is useful to work out the rotation operators around 
the direct lattice points A and B which contain new sym- 
metries not included in the rotation operator around a 
dual triangular lattice point Rn/3 : 

r>A \ s i2it f (ai — 2a\a^— 2a\— 2a^) 

-"■2^/3 — 

x 

I — a 2 + 1, ax — <x 2 >< ax, a 2 | 

pB _ i2 7 r/((a 2 -l) 2 -2aia 2 +2a 1 ) 

2-7T/3 / j 

X 

I - a 2 + l,ax — a 2 + 1 >< ai,a 2 | (10) 
They act on the two vortex fields ip± as: 

(11) 

It is ease to see the two operators Ai, A 2 in Eqnj9]are 
also invariant under Eqn llll Finally we reach the most 
general quartic term invariant under all the above trans- 
formations in Eqns 171 [TH 

U = 7o(|V;+| 2 + IV'-I 2 ) 2 - 7i(l^+| 2 - \i>-\ 2 ? (12) 

In the square lattice, In Landau gauge, a unitary trans- 
formation to the permutative representation is needed to 
reach EqnfTJ] Here in the gauge chosen in Eqn[31 ip± 
are automatically in the permutative representation. An 
important and subtle point is how to construct boson 
density order parameters to characterize the symmetry 
breaking patterns in the direct lattice in terms of the dual 
vortex fields in the dual lattice. In 5 -, the density order 
parameter was constructed to be the most general gauge 
invariant and bilinear combinations of the vortex fields 
ipl. Then it was evaluated at the direct lattice points, 
links and dual lattice points to represent boson density, 
kinetic energy and amplitude of the ring exchanges re- 
spectively. The dual lattice of a square lattice is still a 
square lattice with the same lattice constant. The link 
points of a square lattice also form a square lattice with 
lattice constant y/2. Putting the direct, dual and link 
lattices together forms a square lattice with lattice con- 
stant 1/2. Although like a square lattice, the honeycomb 
lattice is also a bi-partisan lattice consisting of two in- 
terpenetrating triangular sublattices A and B, its dual 
lattice is a triangular lattice which is a frustrated one ( 
Fig.l), its link points form a Kagome lattice. Putting the 
three lattices together forms a very complicated lattice. 
So the density wave order parameters in the honeycomb 
lattice may not be evaluated in the three lattices. Indeed, 
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when trying to identify the density operators in the insu- 
lating state, we find the density operators proposed ir£ 
does not apply anymore in the honeycomb lattice. In 
the following, by studying how gauge invariant bilinear 
vortex fields transform under the complete PSG, we can 
identify both the boson density and boson kinetic energy 
operators in the direct lattice. We need evaluate these 
quantities only at dual lattice points to characterize the 
CDW and VBS orders in the direct lattice. Let's look 
at the generalized density order operators which charac- 
terize the symmetry breaking patterns in the insulating 
states. In the low energy limit, the vortex field is: 

= e -^+- s tp + (x) + e-^-'^_(x) (13) 

Intuitively, the generalized density operator p(x) = 
tp^ (x)ip(x) can be written as: 

p(x) = ^ V+ + ipLiff- + e^iplip- + e~ l & s i>Uj+ (14) 

where Q = 2tt/3(1,1). 

When taking continuum limit on the dual lattice, we 
take one site per dual unit cell which contains two direct 
lattice sites A and B ( Fig. 1). So EqnQj] should contain 
both the information on the boson densities on sites A 
and B and the boson kinetic energy on the link between 
A and B ( or XY exchange energy < S^Sg + h.c. > 
in the spin language ). The single boson Green func- 
tion hopping on the direct lattice is related to the gauge- 
invariant single vortex Green function on the dual lattice 
in a highly non-local way^l. Fortunately, the two boson 
quantities such as the density and the kinetic energy may 
have simple local expressions in terms of the dual vortex 
fields. By studying how the operators in Eqn. [T?] trans- 
form under [71 and I 111 in the scaling limit, we can identify 
these quantities as ( up to a unknown prefactor)22 : 

p A = i>+ip+, pB=ip-ip- 
K AB = e l ^ 3 ^_+e-^ 3 ^+ (15) 

where x stands for dual lattice points only. 

Moving slightly away from half filling / = 1/2 cor- 
responds to adding a small mean dual magnetic field 
H ~ Sf = f — 1/2 in the action. It can be shown that 
inside the SF phase, the most general action invariant 
under all the MSG transformations upto quartic terms is 

a—a/b 

+ ^(e^A^A A -27T(5/V) 2 + U (16) 

where is a non-compact U(l) gauge field. Upto 
the quartic level, with correspondingly defined Vv/b i n a 
square lattice, Eqn[T75]is the same as that in the square 
lattice derived ir£. Because Eqn|l6] is a long wavelength 
effective action, the relations between the phenomenolog- 
ical parameters in Eqn ll6l and the microscopic parame- 
ters in Eqn[T]are not known. Fortunately, we are still able 
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FIG. 1: (a) Bosons at filling factor / are hopping on a hon- 
eycomb lattice ( solid line ) which has two sublattices A and 
B. Vortices are hopping on its dual lattice which is a triangu- 
lar lattice ( dashed line ) which has three sublattices X, Y, Z. 
In the easy-plane limit shown in Fig. 2c, one of the three VBS 
states^ on the honeycomb lattice is shown by the thick bonds 
in the figure. The other two VBS can be obtained by 

or 

to classify some phases and phase transitions and make 
some concrete predictions from Eqn llGI without knowing 
these relations. In the following, we assume r < in 
EqnJSl so the system is in insulating states. We will dis- 
cuss the Ising limit and the easy-plane limit respectively. 

III. ISING LIMIT 

If 71 > 0, the system is in the Ising limit, the mean 
field solution is ip a = 1, ipb = or vice versa. The system 
is in the CDW order. We can see how pa,Pb transform 
under Eqns. 171 and [Til 

R 2tt/3 : PA/B -> PA/B] Rir/3-Ja ■ PA PB (17) 

These transformations confirm that pa and p B indeed 
can be identified as the boson density operators at direct 
sublattices A and B. In this section, we discuss at and 
away from half filling respectively. 

A. SF to CDW transition at half filling Sf = 0. 

If r < 0, the system is in the CDW order which could 
take checkboard (ir, ir) order or a stripe order- 17 . EqnfTol 
is an expansion around the uniform saddle point < V x 
A >= f = 1/2 which holds in the SF and the VBS ( to 
be discussed in section 4 ), so it breaks down in the CDW 
side. So it can not be used to study the SF to the CDW 
transition. 

In the CDW state where r < 0, a different saddle point 
where < V x A a >= 1 — a for sublattice A and < V x 
A b >— a for sublattice B with a < 1/2 should be used. 
It is easy to see that there is only one vortex minimum 
-0b in the staggered dual magnetic field at Sf = 0, the 
effective action inside the CDW state is: 

Cc-cdw = \{d ll -iAl)^ b \ 2 + f\^ b \ 2 + u\^ b \ A + --- 

+ -^{e^xdvAif (18) 
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where f < 0, so the system is in the CDW state where 

< ip b >t^ 0. Eqn[TH]is essentially the same as Eqn lA2l in 
q = 1 case. Note that because < tjj a >^ 0, the gauge field 
A a is always massive, so it does not appear in EanfTSl 
As explained in the appendix, this indicates the density 
fluctuations in sublattice A is suppressed, so can be taken 
as fixed. 

Due to the change of the saddle point structures, the 
transition from the SF to the CDW driven by the horizon- 
tal axis ( quantum fluctuation r ) in Fig. 2b is a strong first 
order transition. If we assumed that Eqn[16] at Sf = 
in the Ising limit could be used to describe this first or- 
der transition, then the CDW side could be described by 

< tpa 0, < Vfc >= or vice versa^S, but no phases 
with < ip a >^ 0, < Vfc 0. However, as shown in 
Eqnmi in the CDW side, < ip a >^ 0, < tp b >^ 0, so 
both gauge fields A a and A b are massive, the density fluc- 
tuations in both sublattice A and B are suppressed, so 
both can be taken as fixed. This fact is due to the change 
of saddle point structure across the SF to the CDW tran- 
sition. So we conclude that due to this change of saddle 
point structure of the dual gauge fields across the SF to 
the CDW, EqnfTn] at Sf — in the Ising limit can not 
be used to describe this first order transition. However, 
it does give a qualitative indication of this strong first 
order transition. This is expected, because in a strong 
first order transition, two separate dual actions Eqnll6l 
in the SF side and Eqn[TH] are needed to describe the two 
sides separately. Note that in the direct picture, the SF 
breaks the U(l) symmetry, the CDW breaks the lattice 
symmetry, so the two sides break two completely different 
symmetries, it can only be first order. 

B. Charge Density Wave (CDW) supersolid away 
from half filling Sf ^ 0. 

Now we look at the effects of the in-commensurability 
Sf = f - 1/2 in Eqn[H In the SF side where r > 0, it is 
known the SF is stable against the change of the chemical 
potential ( or adding bosons ). In the CDW side where 
f < 0, moving slightly away from half filling / = 1/2 
corresponds to adding a small mean dual magnetic field 
H ~ Sf = f - 1/2 in Eqn[H] 

Cic-cdw = 1(8^ - tAl)^ b \ 2 + fy Jb \ 2 + u\^ b \ 4 + ■ ■ ■ 

+ -^(e^xd„A b x -2TTSf5^ T ) 2 (19) 

where the vortices in the phase winding of ip b should be 
interpreted as the the boson number—. Note that be- 
cause < ip a >t^ and remains non-zero in the presence 
of Sf, the gauge field A a is always massive and remains 
massive in the presence of Sf, so it still does not appear 
in EqnrTS] This indicates the density fluctuations in sub- 
lattice A remains suppressed, so can be taken as fixed in 
the presence of Sf. 

EqnrT5]is essentially the same as Eqn lA4l which has the 
structure identical to the conventional q = 1 component 



Ginzburg-Landau model for type-II " superconductors " 
in a "magnetic" field. It was shown in 1 ^ that for type 
II superconductors, the gauge field fluctuations will ren- 
der the vortex fluid phase intruding at H c \ between the 
Messiner and the mixed phase ( see Fig. 2a). For param- 
eters appropriate to the cuprate superconductors, this 
intrusion occurs over too narrow an interval of H to be 
observed in experiments. In the present boson problem 
with the nearest neighbor interaction V± > in EqnQ] 
which stabilizes the (ir,ir) CDW state at / = 1/2, along 
the dashed line driven by the vertical axis ( chemical 
potential p ) in Fig. 2b, this corresponds to a CDW su- 
persolid (CDW-SS) state intruding between the commen- 
surate CDW state at / = 1/2 and the in-commensurate 
CDW state at 1/2 + Sf which could be stabilized by fur- 
ther neighbor interactions in Eqn[T] The first transition 
is in the z — 2,v — 1/2,7] — universality class first 
discussed inA, while the second is a 1st order transition. 
We expect the intruding window at our q = 2 system 
is much wider than that of the q — 1 system. In the 
CDW states, < ip b >^ 0, so the dual gauge field A°^ in 
Eqn[TH] is massive. In the CDW-SS state, < ip b >= 0, 
there is the gapless supcrfluid mode represented by the 
dual gauge field A^. The CDW-SS has the same (tt, tt) 
diagonal order as the C-CDW. The identified universality 
class of the CDW to the CDW-SS transition has many 
physical implications. For example, near the C-CDW 
to the CDW-SS transition, the superfluid density should 
scale as p s ~ \p - i/2\( d + z - 2 > = \p - 1/2| = Sf with 
a logarithmic correction. The logarithmic correction will 
be calculated in a future publication . There must be a 
transition from the CDW-SS to the SF inside the win- 
dow driven by the quantum fluctuation r in the Fig. 2b. 
The universality class of this transition is likely to be first 
order and will be investigated further in a future publi- 
cation . Combining the results in III-A and III-B loads 
to the global phase diagram Fig. 2b in the Ising limit 24 . 



IV. EASY-PLANE LIMIT. 

If 7x < 0, the system is in the easy-plane limit, 
the mean field solution is ip a = e l6a , ip b = e l9b , then 
PA = Pb = 1, so the two sublattices remain equivalent. 
This limit could be reached by possible ring exchange in- 
teractions in Eqn[T] In the following, we discuss at and 
away from half filling respectively. 



A. SF to VBS transition at half filling Sf = 0. 

Because the two sublattices remain equivalent, the uni- 
form saddle point <Vxi>=/=l/2 holds in both the 
SF and the VBS. The system has a VBS order with the 
kinetic energy order parameter Kab = cos(Q ■ x + 8 J) 
where 6- = 9 a — 9 b . Let's look at how the kinetic energy 
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(a) 



FIG. 2: (a) The phase diagram slightly away from / = 1 
described by Eqn lA4l in the appendix is the same as type-II 
superconductors in an external magnetic field H where there 
should be a vortex liquid state intruding between the Messiner 
state and the vortex lattice state. But the intruding regime 
is too narrow to be seen in type-II superconductors^, (b) 
and (c) are the zero temperature global phase diagrams of 
the chemical potential fi versus r in Eqn ll6l in the honeycomb 
lattice, (b) The Ising limit 71 > discussed in section III. 
There is a CDW supersolid (CDW-SS) state intruding be- 
tween the commensurate CDW ( C-CDW ) state at / = 1/2 
and the in-commensurate CDW (IC-CDW) state at 1/2 + 5 f. 
The CDW-SS has the same lattice symmetry breaking as the 
C-CDW. (c) The Easy-Plane limit 71 < discussed in sec- 
tion IV. There is a Valence Bond Supersolid (VB-SS) state 
intruding between the commensurate VBS (C-VBS) state at 
/ = 1/2 and the in-commensurate VBS ( IC-VBS ) state at 
1/2 + 8 f . The VB-SS has the same lattice symmetry breaking 
as the C-VBS. The thin ( thick ) line is the 2nd ( 1st ) order 
transition. The 1st order transition in the Ising ( Easy-plane 
) limit is strong ( weak ) one. (b) and (c) are drawn only 
above half filling, fj, — > — fj, corresponds to below half filling. 



triangular lattice can be divided into three sublattices 
X,Y,Z where x\ + x% — 0,1,2 mod 3, hence Q ■ x — 
0, 27r/3, — 27r/3 respectively ( Fig.l). If we evaluate K at 
X, Y, Z, we find K x = cos 0, K Y = cos(6> + 2tt/3) , K z = 
cos(# — 2ir/3). If we know one bond in one dual unit cell, 
say, labeled by X, then by the translation, we can get 
all the other bonds in sublattices Y and Z with the same 
orientation. The two bonds on the other two orientations 
can be reached by rotations listed in Eqn. [2(3 By this 
way, we can get all the bonds in the whole direct lattice. 
So the system is in the Valence Bond Solid (VBS ) state, 
one VBS is shown in Fig.l. 

If A < 0, 6 = 2mr/3 = 0, 2tt/3, -2tt/3. The kinetic en- 
ergy K = 2,-1,-1 takes also only two values, one strong 
bond and two weak bonds, their ratio is also 2. This case 
is essentially the same as A > 0. Because the sign of the 
kinetic term can be changed in a bi-partisian lattice by 
changing the sign of bi in one of the two sublattices in 
EqnQ] (or can be changed by changing the sign of Sf, Sf 
in one of the two sublattices, but keeping Sf untouched in 
Eqnj2] ) , but the product of the sign around a hexagon is 
fixed. The two cases have the same sign product, so can 
be transformed to each other by the transformation. This 
can also be understood by observing that C3 in Eqnl2"Tl 
behaves like a hopping term in 3 power, so its sign can be 
changed by the transformation. This is in sharp contrast 
to the square lattice to be discussed in section V where 
one sign leads to a columnar dimer, the other leads to a 
plaquette pattern. 



K transform: 

T a : K -> cos(Q • (x - a a ) + 9) 
R„ /3 : K -> cos(Q -x-9); I a : K -> K 

R il% '■ K ~^ cos (Q ■ x + 9 + 27r/3) (20) 

Note that K transforms differently under R w /s and I a , 
because the latter are anti-unitary operators. These 
transformations confirm that K indeed can be identified 
as the boson kinetic energy or the XY exchange energy 
operators at the direct lattice. 

Upto the quartic order, the relative phase between ip + 
and ip- is undetermined. Higher order terms are needed 
to determine the relative phase. It is clear to see there 
are only 3 sixth order invariants: 

Ci = |V + | 6 + |<A-| 6 

c 2 = (|V;+| 2 + |V-l 2 )l^ + l 2 h/>-l 2 

C 3 = A[(^-) 3 + (V>1^+) 3 ] = Acos3# (21) 

where 9 = 9 + — 9-. Especially C3 is invariant under 
R%1% Obviously, only the last term C3 can fix the rel- 
ative phase. This term corresponds to the 3-monopole 
operator in the spin language. 

If A > 0, 9 = (2n + l)7r/3 = tt/3, tt, -tt/3. The kinetic 
energy K = 1,-2,1 takes only two values, one strong 
bond and two weak bonds, their ratio is 2. The dual 



B. Valence Bond Supersolid away from half filling 

Now we look at the effects of the in-commensurability 
Sf = f — 1/2 in Eqnfinj Again, the SF side where r > 
is stable against the change of the chemical potential. In 
the VBS side where r < 0, slightly away from the half- 
filling, EqnfTBl becomes: 

Cvbs = {\d»9 + -A tl f + -^ I (e^xd u A x -2n5f5^y + --- 

+ i\d^-f + 2Acos30_ (22) 

where 6± = 9 a ±9b. 

Obviously, the 6- sector is massive ( namely, 9 a and 9b 
are locked together ) and can be integrated out. Assum- 
ing A > 0, then 6_ = tt. Setting ip + — e l9 + ~ ip a ~ —tpj, 
in EqnfTB] leads to Eqn[l9]with u = 270, so the discus- 
sions on Ising limit case following Eqn llDI also apply. In 
the present boson problem with possible ring exchange 
interactions in Eqn[T] which stabilizes the VBS state at 
/ = 1/2, along the dashed line driven by the vertical 
axis ( chemical potential fi ) in Fig. 2c, this corresponds 
to a VBS supersolid (VB-SS) state intruding between 
the commensurate VBS ( C-VBS) state at / = 1/2 and 
the in-commensurate VBS (IC-VBS) state at 1/2 + Sf 
as shown in Fig. 2c. In the C-VBS state described in the 
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subsection A, < ip a >= — < ipi, >^ 0, < ipli/j a >=< 
V'bV'b >= — < V'lV'b 0, the gauge field is mas- 
sive due to the Higgs mechanism. In the VB-SS state, 

< > = < A >= 0, but < Vl^a > = < i>ti>b >= - < 

iplipb 0, so there is a VBS order characterized by the 
order parameter Kab — cos(Q-x + 9-) which is the same 
as the C-VBS, while the gauge field is massless which 
stands for the gapless superfluid mode inside the VB-SS. 
So this VB-SS has both the VBS order and the super- 
fluid order which justifies its name. So in this VB-SS, the 
density fluctuation on each site is very large signalizing 
its superfluid nature, while the VBS order is fixed signal- 
izing its VBS nature . In this IC-VBS state, Sf valence 
bonds shown in Fig.l is slightly stronger than the others, 
these Sf slightly stronger bonds also form a dilute lattice 
on top of the underlying C-VBS lattice. Again, the first 
transition from the C-VBS to the VB-SS is in the z = 
2, v = 1/2, T] = universality class, while the second from 
the VB-SS to the IC-VBS is 1st order. Near the C-VBS 
to the VB-SS transition, the superfluid density should 



scale as p s ~ \p— 1/2| 



(d+z-2)v _ 



p-l/2\=Sf with log- 



arithmic corrections. The nature of the transition from 
the VB-SS to the SF where < ^ a >=< tp b >= and 
< ip^ipb >= inside the window driven by the quantum 
fluctuation r in the Fig. 2c is likely to be weakly first order 
and will be studied further in a future publication. Com- 
bining (4a) and (4b) leads to the global phase diagram 
Fig.2c. 



V. SQUARE LATTICE. 

As said below EqnfTHl with correspondingly defined 
tpa/b m a square lattice, upto the quartic level, Eqnll6l is 
the same as that in the square lattice derived in 5 . 



A. Ising limit 

So the phase diagram in the Ising limit Fig. 3b remains 
the same as Fig. 2b. The SF to the C-CDW transition at 
Sf = along the horizontal axis is a strong first order 
one. Away from the half filling, along the dashed line in 
Fig. 3b, the IC-CDW can be stabilized only by very long 
range interactions in EqnQ] if it is not stable, then Fig. 3b 
reduces to Fig. 4b ( which is the Fig. 14 inM). 





r<0, y<o 

(c) 



FIG. 3: Phase diagram in square lattice, (a) is identical to 
2a and presented here just for completeness, (b) is similar 
to Fig.2b on the honeycomb lattice in the Ising limit. In the 
easy plane limit (c), there is a possible 2nd order transition 
between the SF and the C-VBS through a so called deconfined 
quantum critical point. However, this has been disputed in 26 . 



component scalar electrodynamics Eqn ll6l at Sf = 
along the horizontal axis is a second order transition 
through a so called "dc-confincd QCP". If this is in- 
deed the case, then C sq — A cos 4$ is irrelevant, there 
is a possible 2nd order transition between the SF and 
the C-VBS through a so called deconfined quantum crit- 
ical point as shown in Fig. 3c. Away from the half fill- 
ing, along the dashed line in Fig. 3c, the phase transi- 
tions are the same as those in the honeycomb lattice 
Fig. 2c. The transition from the VB-SS to the SF where 
< ipa >—< tpb >= and < ip^i/jb >= inside the win- 
dow driven by the quantum fluctuation r in the Fig. 3c is 
also through line of deconfined quantum critical points. 
Recently, some connections are made between the cor- 
relation functions ( which lead to transport properties 
such as conductivity ) near the deconfined quantum crit- 
ical point in Fig. 3c and those in M = 8 supersymmet- 
ric SU{N) Yang-Mills gauge theory in the large N limit 
at d = 2 49 . These correlation functions were calculated 
through AdS A x S 7 /CFT 3 connection^. The particle- 
vortex self-duality in Fig. 3c and corresponding electro- 
magnetic self-duality in AdS^ x S 7 /CFT3 put strong 
constraints on correlation functions in the two models 
respectively. 

However, recently, by using QMC simulations at much 
larger systems, the authors in 26 pointed out that the 
QMC study in 2 ^ violates the hyperscaling and concluded 
that there is no such deconfined QCP, so the transition 
from SF to the VBS driven by the horizontal axis ( quan- 
tum fluctuation r ) is a weak first order one. If this is 
indeed the case, the phase diagram in the easy-plane limit 
remains the same as in a honeycomb lattice Fig. 2c. 



B. Easy-plane limit 

In the Easy-plane limit, as shown ir£, the lowest order 
term coupling the two phases 8 a /b is C sq — A cos AO, so 
the C3 term in Eqn(2H need to be replaced by C sq . If 
A is positive ( negative), the VBS is Columnar dimer ( 
plaquctte ) pattern^. 

By using duality analysis and QMC simulation, the au- 
thor ir* 2 ^ suggested that in the easy plane limit, the q = 2 



VI. IMPLICATION ON QUANTUM 
MONTE-CARLO (QMC) SIMULATIONS. 

The EBHM Eqn[T] of the hard core bosons on square 
lattice with V\ and V% interactions was studied by the 
QMC simulations ir*ii, the authors found a stable striped 
(tt,0) and (0,tt) SS ( Fig.4 ). A stable (7T,7r) SS can 
only be realized in the soft core boson casein. But the 
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M (b) 



FIG. 4: Phase diagrams of the EBHM Eqn[T] in the hard 
core limit U = oo achieved from the QMC in^: (a) with 
the nearest neighbor interaction Vi > 0, the CDW ordering 
wavevector is (ir,n). The corresponding (n,n) supersolid is 
unstable against the phase separation. However, it becomes 
stable in the soft core limil^ (b) Vi = 0, but next nearest 
neighbor interaction V2 > 0, the CDW ordering wavevector is 
(-7T, 0) or (0, tt) which is called stripe phase. The corresponding 
stripe supersolid is stable even in the hard core limit. So there 
is a narrow window of stripe SS intervening between the stripe 
solid and the SF. The universality class of the stripe solid to 
the stripe SS in (b) was not studied in^. The thin (thick) 
line is a 2nd (1st ) order transition. 



nature of the CDW to supersolid transition has never 
been addressed. Our results in the Fig. 3 show that the 
CDW to the SS transition must be in the same univer- 
sality class of Mott to superfluid transition with exact 
exponents z — 2, v = 1/2, 77 = with a logarithmic 
correction. It is important to (1) confirm this predic- 
tion by finite size scaling through the QMC simulations 
in square lattice for (0, 7r) and (it, 0) supersolid in both 
hard core and soft core case in the Fig. 4 and (tt, tt) su- 
persolid in the soft core case only (2) do similar things 
in honeycomb lattice to confirm Fig.l. Some of our re- 
sults in Fig. 3 are indeed confirmed in very recent QMC 
simulations on soft core bosons in honeycomb lattice 2 ^. 
(3) To Eqn[T] with U = 00, V% > 0, adding a ring ex- 
change term -K s J2ijki( b l b j b l b i + /i.e.) where i,j,k,l 
label 4 corners of a square in the square lattice and 
-K h J^tjkimn ( b l b j b l b i b L b n + h.c.) where i,j,k,l,m,n la- 
bel 6 corners of a hexagon in the honeycomb lattice to 
stabilize the C-VBS state at half filling ( Note that if 
K s ,Kh > 0, the QMC are free of sign problems), then 
confirm the prediction on C-VBS to VB-SS transition 
in Fig. 2c and Fig. 3c. We expect that the VB-SS phase 
should be stable against a phase separation in both hard 
and soft core cases. The second transition ( CDW-SS to 
IC-CDW in Fig.2b and Fig.3b and the VB-SS to IC-VBS 
in Fig. 2c and Fig. 3c ) is hard to be tested in QMC, 
because some very long range interactions are needed to 
stabilize the IC-CDW or the IC-VBS state. They are 
first order transition anyway. 

In fact, one of the predictions in this paper on the scal- 
ing of the superfluid density p s ~ \p — 1/2| was already 
found in the striped (tt, 0) solid to striped supersolid tran- 
sition by QMC in Sec.V-B of RefJ^. As shown in section 
III-B, there should be a logarithmic correction to the 
scaling of p s , it remains a challenge to detect the loga- 



rithmic correction by a high precision QMC. Of course, 
the superfluid density is anisotropic p x s > p v s in the (tt, 0) 
supersolid, but they scale in the same way with different 
coefficients 14 . Although the authors ini 4 . suggested it is 
a 2nd order transition, they did not address the univer- 
sality class of the transition. 

In the following three sections, I will discuss the ap- 
plications of the results achieved in the previous sections 
on 3 different experimental systems. 



VII. ULTRA-COLD ATOM ON SQUARE AND 
HONEYCOMB OPTICAL LATTICES 

Recently, Bose-Einstein condensation (BEC) was re- 
alized in ultra-cold atomic gases ( for a review , see^ 8 - ). 
Superfluid to Mott insulator transition was also observed 
in optical lattices of ultra-cold alkali atoms2£. Atomic 
physicists are constructing effective various kinds of 2d 
and 3d optical lattices using laser beams and then load 
either ultra-cold fermion or boson atoms at different fill- 
ing factors on the lattices. They may tune the param- 
eters in Eqn[T] to realize different phases and quantum 
phase transition s 2 ^ 28 . The optical honeycomb lattice ge- 
ometry Fig.l and square lattice could be realized in fu- 
ture ultra-cold atomic experiments. The challenge is to 
achieve longer range interactions than the onsite interac- 
tion. Very exciting perspectives to achieve longer range 
interactions have been opened by recent experiments^ on 
cooling and trapping of polar molecules. Being electri- 
cally or magnetically polarized, polar molecules interact 
with each other via long-rang anisotropic dipole-dipole 
interactions. Loading the polar molecules on a 2d optical 
lattice^ with the dipole moments perpendicular to the 
trapping plane can be mapped to EqnQ]with long-range 
repulsive interactions ~ p 2 /r 3 where p is the dipole mo- 
ment. There are also other ways to generate long-range 
interactions^! We expect Fig. 2a can be easily realized. 
There could also be some efficient ways to generate the 
ring exchange interactions^! needed to stabilize the the 
VBS state in Fig.2b. 



VIII. ADATOM ADSORPTIONS ON 
SUBSTRATES ( HONEYCOMB LATTICE ) 

In this section, we will use the phase diagram Fig. 2a to 
discuss two experimental systems: (1) the reentrant "su- 
perfluid" detected in a narrow region of coverages in the 
second layer of He4 ( called 4 ffe/ 4 -ffe/graphite system 
) adsorbed on graphite by a previous torsional oscilla- 
tor experiment ( 38 . (2) the reentrant "fluidlike" state 
of hydrogen adsorbed on Krypton-preplated graphite ( 
i?2/Kr/graphite ) near half filling which was investigated 
in a recent experiment^. 

The dual vortex approach is a Magnetic Space Group 
( MSG ) symmetry based approach which can be used 
to classify all the possible phases and phase transitions. 
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However, the question if a particular phase will appear 
or not as a ground state depends on the specific values of 
all the possible parameters in the boson Hubbard model 
in EqnflJ namely, the specific values in the real systems, 
it can not be addressed in this approach. There are two 
ways to remedy this short-coming. (1) we can compare 
our theoretical classifications with some known phases 
observed in experiments, then we can be more specific 
on our predictions on the nature of unknown phases and 
phase transitions. (2) As said in section VI, a microscopic 
approach such as Quantum Monte-Carlo (QMC) may be 
needed to supplement the dual field theoretical approach 
. In this section, we will take the first strategy to study 
the boson Hubbard model EqnQ] in honeycomb lattice 
Fig.l near q = 2 ( Fig.l ). The honeycomb lattice in Fig.l 
may describe the preferred adsorption sites in the two 
systems. So the He/ -He/graphite and Z?2/Kr/graphite 
systems may have the same symmetry and belong to the 
same universality class. 

A. 4 He/ 4 He/ graphite system 

A superfluid is a fluid that flows through the tini- 
est channels or cracks without viscosity. So far, the 
phenomenon of superfluidity has been firmly observed 
in only two kinds of systems. The first system is the 
two isotopes of Helium: 4 He and 3 He. He4 is a boson 
and becomes a superfluid when T < T c = 2.18K, while 
He3 is a fermion, two He?> atoms start to pair up when 
T <T C — 2AmK and form a superfluid 3 ^ 3 ^. The second 
system is Bose-Einstein condensation in ultra-cold alkali 
atomic gases ( for a review, see^ ). Recently, intensive 
research activities have been lavished on searching for 
excitonic superfluid in excitons in electron-hole semicon- 
ductor bilayers^i. Torsional oscillator method was used 
to study A He films adsorbed on graphite 38 , a large Non- 
Classical Rotational Inertial ( NCRI ) was detected in 
a narrow window of coverages c = 17 ~ 19 atoms /nm 2 
in the second layer ( 4 iJe/ 4 i?e/graphite system ). The 
NCRI is a low temperature reduction in the rotational 
moment of inertia due to the superfluid component of 
the state. This important phenomenon was interpreted 
as reentrant "superfluid" in this narrow window- 3 ^. Re- 
cently, also by using the torsional oscillator measurement, 
a PSU group lead by Chan observed a marked 1 ~ 2% 
NCRI even in bulk solid 4 He at ~ 0.2K&. This was inter- 
preted as a new state of matter called supersolid which 
has both superfluid order and crystalline order—. How- 
ever, so far, no NCRI in bulk solid H2 was detected. 

At the completion of the first layer, the 4 He atoms 
with the coverage c ~ 12 atoms /nm 2 form a triangu- 
lar lattice which is incommensurate with the underlying 
graphite. It is reasonable to assume the lattice structure 
of the preferred adsorption sites on the second layer form 
a honeycomb lattice ( Fig.l). When applying the EqnQ] 
to the 4 -ffe/ 4 -ffe/graphite system in Fig.l, b\ is the 4 He 
atom creation operator. The hopping amplitude t is de- 



termined by the trapping potential from the incommen- 
surate first 4 He layer. We assume the 4 He atoms interact 
with each other with the LJ potential. Note that the LJ 
potential is quite accurate in describing the long distance 
attractive part, but very crude in the short distance re- 
pulsive part. The repulsive part in real case is expected 
to be softer. The chemical potential fj, determines the 
coverage. The filling factor / is related to the coverage 
c of 4 He with c ~ 19.55 atoms/nm 2 corresponding to 
/ = 1/2 where one 4 He atoms occupies every two lattice 
sites of the honeycomb lattice ( for example, all the A 
sublattice sites ) to form a close packed triangular lattice 
with <Iaa ~ 4.56A and <Iab ~ 2.63A which is just smaller 
than Lennard- Jones (LJ) parameter a( 4 He) ~ 2.64A. 
The onsite U is very big. Because cLab < ^"-He < d-AA, 
Vi is positive, while V2 and further neighbor interactions 
are weakly attractive and can be neglected. Note that 
our theoretical value 19.55 atoms/nm 2 is only 3% away 
from the experimental value 19 atoms/nm 2 . Fig. 2a is 
very similar to the p s versus c ~ / phase diagram in 
4 -ffe/ 4 -ffe/graphite structure near c = 19 atoms/nm 2 ( 
see Fig. 1 in 3 ^ ). The dashed line is the experimental 
path in 4 f/e/ 4 £fe/graphite at T = 0. Ir shows that a 
reentrant supersolid (SS) state is a generic state sand- 
wiched between the C-IC transition. The reentrant " 
superfluid " state in the second layer of the 4 He films 
adsorbed on graphite could be a 4 He supersolid state. 
Very interestingly, the data in the torsional oscillator ex- 
periment in^ 8 - do not show the characteristic form for a 
2d 4 He superfluid film, instead it resemble that in 8 - char- 
acteristic of a possible supersolid in terms of the gradual 
onset temperature of the NCRI, the unusual tempera- 
ture dependence of T$s on the coverage. From Fig. 2a, 
it is easy to see that it maybe difficult to reach the su- 
perfluid by moving along the horizontal r axis, but it is 
very easy to get to the CDW-SS state by moving along 
the vertical coverage axis. In principle, there should be a 
Kosterlitz-Thouless (KT) finite temperature phase tran- 
sition above which the CDW-SS becomes a coexistence 
of the CDW state and normal fluids of interstitials. The 
effects of disorders on this finite temperature transition 
will be studied in and compared with Fig. 2 and 3 ir*2£. 



B. //2/Kr/graphite system 

There are great interests in finding superfluidity in 
other substances. Hydrogen molecules are relatively 
light, quantum fluctuations are large at very low tem- 
perature, a para-H2 molecule is a boson and very similar 
to a He4 atom, in principle, H2 could become superfluid 
at low temperature as He4 docs. Unfortunately, unlike 
He4, due to deeper attractive potentials, bulk H2 solidi- 
fies at low temperature ( T c ~ 1AK ), this preempts the 
possible observation of the speculated superfluity. At suf- 
ficiently high pressure, the solid hydrogen will transform 
into a metallic alkali-like crystal or an unusual two com- 
ponent ( protons and electrons ) quantum liquid at low 
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temperatures? . Potential avenues to prevent p—H2 from 
being solidified to low enough temperature such that the 
superfluid behaviour can be observed is by decreasing the 
average number of neighbors of a H2 molecule. Monte 
Carlo simulations of H2 n clusters with n = 13 — 18 
showed that there is indeed a tendency for the system to 
be superfluid below IK— Indeed, Infra-red spectroscopy 
also provided evidences for superfludity in o(Z?2)n and 
p{H2)n clusters with n — 14 — 17 36 . Another route is 
by the reduction of dimensionality. Extensive theoretical 
and experimental work has been devoted to study H2 and 
D2 films in a variety of substrates. Path Integral Monte- 
Carlo (PIMC) simulations indicated that introduction of 
certain impurities can stabilize a 2 dimensional liquid 
hydrogen which undergoes a Kosterlitze-Thoulcss (KT) 
transition below 1.2 K. 

In a recent experiment, neutron scattering measure- 
ments were used to characterize all the possible phases 
of Z?2 coadsorbed on graphite preplated by a monolayer 
of Kr called f^/Kr/graphite structure 4 ^. Because two 
dimensional graphite has the honeycomb net structure, 
the corresponding preferred adsorption sites form a tri- 
angular lattice. The precoated Kr atoms will occupy the 
preferred adsorption sites and form a triangular lattice. 
Then the D2 deposited on top of the Kr monolayer will sit 
on the preferred adsorption sites of the triangular lattice 
of the Kr monolayer to form a honeycomb lattice. So, the 
lattice geometry is very similar to Fig.l with Kr atoms 
sitting on the triangular lattice, while the H2 molecules 
hopping on the honeycomb lattice. The filling factor / is 
related to the coverage c of D 2 - In the coverage c verse 
the temperature T phase diagram, an unusual feature 
is that in a small coverage range ( 1.20 < c < 1.25 ) 
at the commensurate-incommensurate (C-IC) transition, 
a reentrant fluid phase squeezes in between the C and 
IC phases down to T = 1.5 K which is the lowest tem- 
perature a liquidlike phase of D2 has ever been found. 
The filling factor / is related to the coverage c of H2 
with c = 1.2 corresponding to / = 1/2 where one if 2 
molecule occupies every two lattice sites of the honey- 
comb lattice with &aa ~ 3.87 A and dAB ~ 2. 24 A in 
Fig.l. Because dAB < &h 2 < 3aa, V\ could also be 
positive and very big, while V2 and further neighbor in- 
teractions are weakly attractive, the above discussions on 
4 ife in subsection A can also be applied to this system. 
Fig. 2a is very similar to the c ~ / versus low T phase 
diagram in f^/Kr/graphite structure near c = 1.20 ( 
Fig. 6 in 4 ^ ). The dashed line is the experimental path 
in ff 2 /Kr/graphite at T = 0. The so called (1 x 1)[±] 
commensurate phase in^S is the CDW state where the 
D2 atoms occupy one of the two sublattices of the un- 
derlying honeycomb lattice. The transition at zero tem- 
perature is a C-CDW to CDW-SS to IC-CDW transi- 
tion described by Fig. 2 (a). The first transition is in the 
z = 2, v = 1/2,1] = universality class. Because v < 1, 
from Harris criterion, disorders are relevant. The sec- 
ond transition is first order. It suggests the reentrant " 
fluidlike " phase could be the H2 supersolid. If this is 



indeed the case, this may lead to the first observation 
of if2 supersolid. Torsional oscillator experiment simi- 
lar to that used inM may be used to measure the NCRI 
of this ff 2 supersolid state. Obviously, this experiment 
can avoid the difficulty of using D 2 which has large co- 
herent neutron scattering cross section, but smaller de 
Boer quantum number. In fact, Kr may not be the best 
spacer to observe the H2 supersolid state. The ideal sit- 
uation is to make dAB is just slightly smaller than an 2 ■ 
Larger atoms like Xe, CF4, SF 6 may be more suitable 
spacers, because they may increase the distance dAB hi 
the Fig.l just slightly smaller than <jh 2 . We suggest 
that 4 ife supersolid should also appear in, for example, 
4 ife/Kr/graphite structure. 

The reentrant lattice SS discussed in this section is 
different from the bulk 4 ife SS state discussed in£, al- 
though both kinds of supersolids share many interesting 
common properties. In the former, there is a periodic 
substrate or spacer potential which breaks translational 
symmetries at the very beginning. While in the latter, 
the lattice results from a spontaneous translational sym- 
metry breaking, so the theory developed in this paper on 
lattices is completely different from the continuum the- 
ory developed ins. Combined with the results ir*2, we 
conclude that 4 ife supersolid can exist both in bulk and 
on substrate, while although f?2 supersolid may not exist 
in the bulk, but it may exist on wisely chosen substrates. 
Ultra-cold atoms supersolid could also be realized in op- 
tical lattices. 



IX. COOPER-PAIR SOLID, SUPERFLUID AND 
COOPER-PAIR SUPERSOLID IN HIGH 
TEMPERATURE SUPERCONDUCTORS ( 
SQUARE LATTICE ) 

Very recently 4 ^, by using both angle resolved 
photoemission (ARPES) and scanning tunneling mi- 
croscopy (STM) on high temperature superconductor 
La2- x Ba x Cu0 4 42 i 43 , Valla et.al detected a quasi-particle 
energy gap with d-wave symmetry even when the super- 
conductivity is completely suppressed at x = 1/8, while 
having almost equally strong superconducting phases at 
both higher and lower dopings ( Fig. 5). The gap turns 
out to reach maximum at x = 1/8 despite T c — ► at 
x = 1/8. This fact suggests that the most strongly bound 
Cooper pairs at x — 1/8 are most susceptible to the 
charge density wave ( CDW ) ordering which suppresses 
T c — > 0. In this section, I propose that the formation of a 
stripe Cooper pair supersolid may be able to explain the 
very unusual phase diagram of La,2- x Ba x CuO± at and 
near doping level x = 1/8. 

The fact that the STM and ARPES measurements 
in 4 ^ still detected an energy gap at Fermi surface with 
d x 2_ y 2 wave symmetry suggests that there are tightly 
bound Cooper pairs in this cuprate. We can treat these 
Cooper pairs as bosons, therefore many important re- 
sults achieved in previous sections on interacting bosons 
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on its own and is also the starting point to incorporate 
quasi-particle a 45 ' 46 ! 47 and spin excitations 4 ^ into the pic- 
ture. As shown in section IV, a valence-bond supersolid 
( VB-SS ) can be stabilized if there is a considerable ring 
exchange interaction. An interesting question to address 
is if valence-bond Cooper-pair supersolid can be realized 
in some of these cuprates. 



o.l 1/8 0.2 o.3 x 

FIG. 5: The dashed line is the energy gap Ao which reaches 
minimum at x — 1/8, the solid line is the superconducting 
transition temperature T c which has two domes and is zero 
at x = 1/8. 



hopping on lattices maybe applied to this cuprate. One 
important new feature for tightly bound Cooper pairs is 
that they carry charges 2e, therefore interact with bare 
long-range Coulomb interaction ( in high T c cuprates, by 
taking the penetration depth A —> oo, one can neglect 
the very weak Meissner effect )£. I propose the following 
physical picture: at exactly x = 1/8, the filling factors 
of the bosons on the square lattice is / = ijp = 7/16, 
the system is a CDW ( more specifically, stripe ) of Cop- 
per pairs with 4 sites per unit cell, so T c = 0. Then we 
need q = 16 dual vortices to describe the stripe CDW in 
the dual vortex picture. It was explicitly pointed out in 
section III that special care is needed to choose the cor- 
rect saddle point of the dual gauge field to describe 
this stripe CDW state. Slightly away from x = 1/8, 
on both sides, the groundstate is a Stripe Cooper pair 
supersolid ( CP-SS ) ! When x > 1/8 (x < 1/8), 
it is a hole ( electron ) doped stripe CP-SS. Indeed, 
as shown by Quantum Monte-Carlo simulation 1 - 4 , only 
striped supersolid is stable in a square lattice for hard- 
core bosons. The quantum phase transition from the 
stripe Cooper pair solid to the stripe CP-SS driven by 
the doping is in the same universality class as that from 
a Mott insulator to a superfluid, therefore have exact 
exponents z = 2, v = 1/2, rj = ( with possible logarith- 
mic corrections ). The resulting stripe CP-SS has the 
same lattice symmetry breaking patterns as the stripe 
Cooper pair solid with the superfluid density scaling as 
p s ~ \x - \/%\(d+z-2)v = ja; _ 1/8| = \5x\ with a possible 
logarithmic correction. As pointed out in section III, the 
superfluid density is anisotropic with p s along the stripe 
is larger than that normal to the stripe, but both scale in 
the same way with different coefficients. From Uemura 
relation, we conclude that T c ~ \5x\ ( in fact, it scales as 
the smaller p s in the Stripe CP-SS ). This result very nat- 
urally explains why T c indeed looks linear in \8x\ and the 
phase diagram is nearly symmetric near x = 1/8 found 
mil. 

The above picture only involves the sector of the 
tightly bound Copper pair. Of course, understanding the 
ground state of the Cooper-pair sector is very important 



X. CONCLUSIONS 

By using the DVM, we studied superfluid, solid and su- 
persolid and quantum phase transitions of the extended 
boson Hubbard model near half filling on bipartite op- 
tical lattices such as honeycomb and square lattice and 
mapped out the global phase diagram at T = in a 
unified scheme. We identified boson density and boson 
kinetic energy operators in terms of the dual vortex fields 
to characterize symmetry breaking patterns in the insu- 
lating states and supersolid states. In the DVM, starting 
from the featureless superfluid state where the the aver- 
age value of the dual gauge field through a dual plaquette 
is taken to be uniform and equal to the boson density 
/ = p/q, its fluctuation is coupled to q dual vortex or- 
der parameters and is gapless, one study all the possible 
symmetry breaking patterns by condensing the q vortex 
order parameters. We first study the transition driven by 
the ratio of the kinetic energy over the potential energy 
at the commensurate fillings / = p/q along the horizonal 
axis in Fig. 2-4. In the Ising limit, we found that the 
saddle point of the dual gauge fields should be chosen in 
a self consistent way in the CDW side, this is in sharp 
contrast to the more familiar / = n case where the self- 
consistency condition is automatically guaranteed. Due 
to this change of saddle point structure of the dual gauge 
fields on the both sides, the SF to the CDW transition is 
always a strong first order one. In the Easy-plane limit, 
we found that the self-consistency condition is automati- 
cally guaranteed, the SF to the VBS transition is a weak 
first order one. Then we study the transition driven by 
chemical potential slight away from the commensurate 
fillings / = p/q along the vertical axis in Fig. 2-4. Wc 
found that in the insulating side, the transition at zero 
temperature driven by the chemical potential must be a 
C-CDW ( or C-VBS ) at half filling to a narrow window 
of CDW- ( VB-) supersolid, then to a IC-CDW ( IC- 
VBS ) transition in the Ising ( easy-plane ) limit. The 
valence bond supersolid is a novel kind of supersolid first 
proposed in this paper. Although the density fluctuation 
at any site is very large indicating its superfluid nature, 
the boson kinetic energies on bonds between two sites 
are given and break the lattice translational symmetry 
indicating its valence bound nature. The first transition 
is in the same universality class as that from a Mott in- 
sulator to a superfluid driven by a chemical potential, 
therefore have exact exponents z = 2,v — 1/2, rj = 
with a logarithmic correction. The second is a 1st order 
transition. The results achieved in this letter could guide 
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QMC simulations to search for all these phases and con- 
firm the universality class of the transitions. This VB-SS 
should be stable in both hard core and sift core limit if 
there is a sufficient strong ring exchange term in EqnfT] 
It will be extremely interesting to search this novel kind 
of supersolid in a specific microscopic model by QMC 
simulation. 

In the second part of this paper, we study the appli- 
cations of the results achieved in the first part by the 
DVM to 3 very important ultra-cold atomic and con- 
densed matter experimental systems. The EBHM in 
honeycomb and square lattices could be realized in ul- 
tracold atoms loaded on optical lattices. So the results 
achieved in this paper may have direct impacts on the 
atomic experiments in optical lattices. Then we applied 
the results to two condensed matter experimental sys- 
tems: (1) adatom adsorption on different substrates such 
as the adsorption in the second layer of 4 He adsorbed on 
graphite and Hydrogen adsorbed on Krypton-preplated 
graphite in honeycomb lattice, we find that a reentrant 
He supersolid in the second layer on graphite may be 
responsible for the NCRI detected in the torsional oscil- 
lator experiments ir*2£ and a reentrant H2 supersolid at 
T = maybe responsible for the reentrant liquid state 
squeezed between C and IC phases detected by coherent 
neutron scattering. We propose that a judicious choice 
of substrate may also lead to an occurrence of hydro- 
gen lattice supersolid. (2) superconducting phase dia- 
gram near x = 1/8 in high temperature superconduc- 
tor La,2-xBa x CuO \ in square lattice. We conclude that 
there is a stripe CDW at 1/8 which suppresses T c to 
zero. There are hole and electron doped Cooper pair su- 
persolids on both sides of x = 1/8 which has the same 
lattice symmetry breaking pattern as the stripe CDW at 
x = 1/8 and with p s ~ T c ~ \x - 1/8|. Although the 
SS in lattice models is different from that in a contin- 
uous systems, the results achieved in this paper on lat- 
tice supersolids may still shed some lights on the possible 
microscopic mechanism and phenomenological Ginsburg- 
Landau theory of the possible 4 He supersolids^. 
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No. PHY99-07949 at KITP-C by the Project of Knowl- 
edge Innovation Program (PKIP) of Chinese Academy of 
Sciences. I also thank CP. Sun for hospitality during my 
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APPENDIX A: DUALITY AT INTEGER 
FILLINGS / = n, THE ROLE OF THE DUAL 
GAUGE FIELD AND THE SELF-CONSISTENCY 
CONDITION 

It is instructive to review the direct boson picture and 
the duality transformation to the dual vortex picture in 
this simplest case with integer filling / = n 4 . The con- 



tents of this appendix are not new. However, we stress 
the self-consistence check on the average value of the dual 
magnetic field on both sides of SF and Mott insulator. 
We also explicitly spell out the physical significance of 
the fluctuations of the dual gauge field on both sides. 
These clarifications are very helpful to motivate the self- 
consistence condition on the average value of the dual 
magnetic fields in the translational symmetry breaking 
insulating sides, especially in the CDW side at q > 2 
cases first discovered ini£. Although this self-consistence 
condition is automatically satisfied at / = n case, they 
become a non-trivial constraint on a self-consistent the- 
ory at q > 2 cases. In a recent unpublished note^, the 
author found this self-consistent condition is even more 
non-trivial and important in frustrated lattices such as 
triangular and Kagome lattices than the bipartite lat- 
tices discussed in this paper. As shown in the main text, 
in both the CDW side ( using limit ) and the valence 
bond ( easy plane limit ), at and slightly away from the 
commensurate filling f — p/q with q > 2 cases, the effec- 
tive actions can be mapped to Eqn lA2l and Eqn lA4l 

In the direct picture, it is convenient to start from the 
Mott insulating side and look at its low energy excita- 
tions, it is easy to see that the creation of a particle by 
(jy is always accompanied by a creation of a hole <j> due 
to the particle hole symmetry at / = n, so the Ginzburg- 
Landau theory to describe the Mott to superfluid tran- 
sition in terms of the boson order parameter <p is given 
by the well know 2 + 1 dimensional relativistic complex 
scalar theory: 

S b = J d 2 rdT[\d T <t ) \ 2 + \S/(i > \ 2 +r\4>\ 2 +u\4>\ 4 + ---] (Al) 

In the Mott state r > 0, so (4>) — 0, there is a Mott 
gap. In the superfluid state r < 0, so {(f)) + 1 0. Due to 
the symmetry breaking in the SF state, there is a gapless 
goldstone mode given by the phase fluctuation of <f>. 

In the dual vortex picture, it is convenient to start from 
the Superfluid side and look at its low energy excitations. 
There are two kinds of low energy excitations. The first is 
just the gapless Goldstone mode in the phase of <p which 
is given by a dual gauge field fluctuation . The second 
is the topological vortex excitation in the phase winding 
of ip. Obviously, the number of vortex ip* is equal to the 
number of anti- vortex ip, so the Ginzburg-Landau theory 
to describe the superfluid to the Mott transition in terms 
of the dual vortex order parameter ip is given by the 2 + 1 
dimensional scalar electrodynamics: 

S d = Jd 2 rdr [ \(d^ - tA^\ 2 + +rM 2 + u d \iA 4 + ■ ■ ■ 

+ \(e^xd v A x ) 2 ] (A2) 

In the superfluid state > 0, so (ip) = 0, there is a 
gapless fluctuation given by the dual gauge field A mu . 
Integrating out the gauge field fluctuation A^ will lead 
to a logarithmic interaction between the vortices. In the 
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Mott insulating state < 0, so (ip) 7^ 0. Due to the 
"symmetry" breaking in the Mott insulating state, the 
dual gauge field acquires a mass due to Higgs mech- 
anism, so there is a Mott gap in the Mott phase. The 
vortex action Eqn lA2l is dual to the boson action Eqn IAll 
Both actions lead to equivalent description of the SF to 
Mott transition which is in the 3DXY universality class 
and the same excitation spectra in both phases. Indeed, 
this universality class has been confirmed by QMC in 
the first reference in 4 - from both the direct boson action 
Eqn IAll and the dual vortex action Eqn lA2l 

In the direct boson picture, a vortex is a singularity 
in the boson wavefunction, so a boson wavefunction ac- 
quires a 2ir phase when in encircles a vortex. In the dual 
vortex picture, a boson is a singularity in the vortex wave- 
function, so a vortex wavefunction acquires a 2n phase 
when in encircles a boson. So the average strength of the 
dual magnetic field of the gauge field A^ through a dual 
plaquette is equal to the boson density f = n, because 
2nn is equivalent to 0, so the average value can be simply 
taken to be zero. It is important to stress that the aver- 
age density of bosons is the same in both the SF and the 
Mott insulating side, namely, it takes the integer n on 
both sides, so the average strength of the dual magnetic 
field can be taken as zero on both side, then the fluctua- 
tions in the dual gauge field reflects the fluctuations of the 
boson density. As said in the last paragraph, the gauge 
field is gapless in the SF phase, so the boson density fluc- 
tuation in the SF is very large, this is expected, because 
the SF is a phase ordered state, so it has a large den- 
sity fluctuation. However, in the Mott phase, the gauge 
filed is massive, so the density fluctuation is suppressed 
in the Mott phase. This is expected also, because the 
Mott phase is a density ordered phase, so has very little 
density fluctuations. So by looking at the behaviors of 



the dual gauge field on both sides, one can distinguish 
the properties of the two phases. In short, although the 
average value of the dual gauge field is the same on both 
the SF and the Mott insulating side, its fluctuations are 
completely different which, in turn, is determined by the 
average value of the vortex order parameter ip in Eqn lA2l 
When slightly away from / = n which is at a in- 
commensurate density, there is no particle-hole symme- 
try anymore, so there should be a first order imaginary 
time derivative, Eqn IAll becomes: 

S lc -b = J d 2 rdT[tfd T (j> + |V</>| 2 + r#| 2 + u\^\ 4 + ■■■} 

where we have dropped the second derivative term |<9 T </>| 2 
which is less important than the linear derivative term. 
If there is only a onsite interaction in Eqn[T] then away 
from / = 1, the system is always in the superfluid state. 
The transition from the Mott to the SF transition driven 
by the chemical potential has the critical exponents z = 
2, v = 1/2, T) = with a logarithmic correction. 

In the dual vortex picture, the linear time derivative 
term corresponds to adding a small mean dual magnetic 
field Sf = f - n to EqnEJ 

S lc - d = J d 2 rdr [ \(d f t-iA^\ 2 + +rM 2 +u d \i>\ 4 + --- 

+ \(e^xd,A x -2ndf5^ T ) 2 ] (A4) 

This action is essentially the same as the Ginzburg- 
Landau model for a superconductor in a external mag- 
netic field if we identify the r direction as the z direction 
along which the magnetic field Sf is applied. For a type 
II superconductor, its phase diagram is shown in Fig. 2a. 
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